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ABSTRACT 

Axisymmetric, orbit-based dynamical models are used to derive dark matter scaling relations for Coma 
early-type galaxies. From faint to bright galaxies halo core-radii and asymptotic circular velocities increase. 
Compared to spirals of the same brightness, the majority of Coma early-types - those with old stellar popula- 
tions - have similar halo core-radii but more than 2 times larger asymptotic halo velocities. The average dark 
matter density inside 2 reft decreases with increasing luminosity and is 6.8 times larger than in disk galaxies of 
the same Z?-band luminosity. Compared at the same stellar mass, dark matter densities in ellipticals are 13.5 
times higher than in spirals. Different baryon concentrations in ellipticals and spirals cannot explain the higher 
dark matter density in ellipticals. Instead, the assembly redshift (1h-z) of Coma early-type halos is likely about 
two times larger than of comparably bright spirals. Assuming that local spirals typically assemble at a redshift 
of one, the majority of bright Coma early-type galaxy halos must have formed around z « 2-3. For about half 
of our Coma galaxies the assembly redshifts match with constraints derived from stellar populations. We find 
dark matter densities and estimated assembly redshifts of our observed Coma galaxies in reasonable agreement 
with recent semi-analytic galaxy formation models. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: formation — galaxies: halos — galaxies: 
kinematics and dynamics — (cosmology:) dark matter 



1. INTRODUCTION 

Present-day elliptic al galaxies are known to host mostly 
old stellar populations (Tr ager et al.l2000llTerlevich & ForbesI 
12002, Thomas et al. 2005a|). Whether their stars have formed 
in situ or whether ellipticals assembled their present-day mor- 
phology only over time (for example by mergers) is less clear 
An important clue on the assembly redshift of a galaxy is pro- 
vided by its dark matter de nsity. For example, in the sim- 
ple spherical collapse model jGunn & Gottll 1972b the average 
density of virialized halos is proportional to the mean den- 
sity of the universe at the formation epoch: halos which form 
earlier become denser. Similarly, in cosmological A^-body 
simulations, the concentration (and, thus, the inner density) 
is found to be higher in halos that have assembled ea rlier (e.g. 
iNavarro. Frenk & WhitdlT996l iWechsler et al.ll2002l) . In ad- 
dition to this connection between formation epoch and halo 
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density, the final halo mass distribution also depends on the 
interplay between dark matter and baryons during the ac- 
tual galaxy formation proces s (e.g. iBlumenthal et alJll986l 
iBinnev. Gerhard & Silkll200ih . Then, the properties of galaxy 
halos provide valuable information about when and how a 
galaxy has assembled its baryons. 

Despite its cosmological relevance, the radial distribution 
of dark (and luminous) mass in early-type galaxies is not 
well known: because of the lack of cold gas as a dynam- 
ical tracer, masses are difficult to determine. Stellar dy- 
namical models require the exploration of a galaxy's or- 
bital structure and have only recently b ecome available for 
axisymmetric or mor e general systems dCretton et al.] 
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dPersic. Salucci & Stel Il996allbl and iKormendv & Freema"ril 
l2004l) . The aim of the present paper is to provide empiri- 
cal scaling relations for generic cluster early-types (flattened, 
with different degrees of rotation). In particular, this paper is 
focussed on the inner dark matter density and its implications 
on the assembly redshift of elliptical galaxy halos. 

The paper is organized as follows. In Sec.|2] we review the 
galaxy sample and its modelling. Dark matter scaling rela- 
tions are presented in Sec. [3] Sec. |4]is dedicated to the dark 
matter density. The effect of baryons on the dark matter den- 
sity is discussed in Sec.|5] while Sec. |6]deals with the halo 
assembly redshift. In Sec.|7]our results are compared to semi- 
analytic galaxy formation models. A summary is given in 
Sec. [8] In the following, we assume that the Coma cluster is 
at a distance of li = lOOMpc. 

2. GALAXY SAMPLE, MODELS AND BASIC DEFINITIONS 

The dark halo parameters discussed in this paper are de- 
rived from the axisymmetric, orbit-base d dynamica l model s 
of bright Coma galaxies presented in Thom as et al.l (l2007ah . 
The original sample comprises two cD galaxies, nine ordinary 
giant ellipticals and six lenticular/intermediate type galaxies 
with luminosities between Mg = -18.79 and Mb = -22.56. 
The spec troscopic and photometric observa t ions are dis- 
cussed in J0rgensen, Franx & Kiaergard ('1996'),'Mehle rt et"an 
(I2OOO), Wegneretal. (2002) and Corsini et al. (2008). Our 
implementation of Schwarzschild's (1979) orbit superposi- 
tion techniqu e for axisymrn etric potentials is described in 
iThomas et al] (|2004 l2005bl) . For a detailed discussion of 
all the galaxy models the reader is referred to iThomas et alJ 
(HoOTi. 

Three of the 17 galaxies from IThomas et akl (l2007al) are 
excluded from the analysis below. Firstly, we do not con- 
sider the two centi-al cD galaxies (GMP2921 and GMP 3329; 
GMP numbers from iGodwin. Metcalfe & PeachI 11983.) . be- 
cause their dark matter profiles may be affected by the cluster 
halo. Secondly, we omit the E/SO galaxy GMP1990, whose 
mass-to-light ratio is constant out to 3 reff. The galaxy either 
has no dark matter within this radius, or its dark matter density 
follows closer the stellar light profile than in any other Coma 
galaxy. In either case, the mass structure of this object is dis- 
tinct from the rest of the sample galaxies. In addition to the 
remaining 14 galaxies we consider two further Coma galaxies 
for which we collected data recently (GMP3414, GMP4822). 
The models of these galaxies are summarized in App.lAl 

Similar dynamical models as used here have been applied to 
the inner regions of ellipticals, where it has been assumed that 
mass follows light (e.g. Gebhardt et alJ2003llCappellari et alJ 
|2006|) . In contrast, our models explicitly include a dark mat- 
ter component (cf. IThomas et al.ll2007a[) . We probed for two 
parametric profiles. Firstly, logarithmic halos 



TABLE 1 

Galaxy parameters shown in Fig.|T] 



PDM(r) = 



4ttG {rj + r^f 



(1) 



which posses a constant-density core of size r/, and have an 
asymptotically constant circular velocity v/,. The central den- 
sity of these halos reads 



Ph = 



Secondly, NFW-profiles 
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galaxy 
(1) 


(2) 


(3) 


1 ''h 

(4) 


(5) 


0144 


10.61 


11.56±0.12 


0.64 ±0.31 


2.33 ± 0.10 


0282 


10.46 


11. 60 ±0.12 


1.23 ±0.24 


2.70±0.10 


0756 


10.89 


11.13±0.12 


1.10 ±0.09 


2.33±0.10 


1176 


10.31 


10.73±0.13 


0.53 ±0.1 8 


2.30 ±0.11 


1750 


10.75 


11 CO 1 rv 1 

11.58ih0.12 


1.2/ ±0.95 


n TA 1 A 

2. /O ± 0.23 


2417 


10.60 


11.43±0.12 


1.38±0.59 


2.70 ±0.38 


2440 


10.30 


11.23±0.12 


1.04±0.15 


2.68±0.13 


3414 


10.13 


11.02±0.15 


0.99 ±0.54 


2.55 ±0.26 


3510 


10.34 


11.28±0.13 


1.07 ±0.39 


2.46 ±0.21 


3792 


10.58 


11.56±0.13 


1.18±0.35 


2.74 ±0.22 


3958 


9.70 


10.81 ±0.13 


0.83 ±0.35 


2.44 ±0.29 


4822 


10.70 


11.69±0.16 


1.11 ±0.50 


2.74 ±0.37 


4928 


11.08 


12.06±0.14 


1.46 ±0.39 


2.71 ±0.19 


5279 


10.72 


11.59±0.12 


1.45 ±0.48 


2.68 ±0.28 


5568 


10.79 


11.89±0.12 


1.82 ±0.39 


2.81 ±0.20 


5975 


10.47 


11.04±0.12 


0.23 ±0.30 


2.30±0.11 


Note. 




- (1) 


Galaxy 


id from 



IGodwin. Metcalfe & PeacU ( Il983h : (2) galaxy B-band 
luminosity Lb; (3) stellar mass M, in solar units; (4,5) 
logarithmic-halo core-radius logr;,/kpc and circular 
velocity log v;, /km s"' . 

which are found in cosmological A^-body simulations 
(|Navarro."F renk & White 1996). The majority of Coma 
galaxies are better fit with logarithmic halos, but the signif- 
icance over NFW halo profiles is marginal. Even if the best fit 
is obtained with an NFW-halo, then the inner region s are still 
dominated by stellar mass (cf. IThomas et al.ll2007ah . In this 
sense, our models maximize the (inner) stellar mass. 

In Sec.|3]we will only discuss results based on logarithmic 
halo fits (i.e. we use the halo parameters from columns 5 and 
6 in Tab. 2 of IThomas et al.l l2b07a and columns 3 and 4 in 
Tab. m of App. IaI respectively). While these are not neces- 
sarily the more realistic profiles, they minimize systematics 
in the comparison with published scaling relations for spirals 
that were performed using cored profiles similar to our loga- 
rithmic halos. The NFW fits are used in Sec.|5] 

The B-band luminosities of Coma galaxies used in this pa- 
per are taken from Hyperleda. We adopt a standard uncer- 
tainty of AMb = 0.3 to account for zero-point uncertainties, 
systematic errors in the sky subtraction, seeing convolution, 
profi le extrapolation and others. Effecti ve r adii are taken 
from 'J 0rgensen. Franx & Kiaergard! (1 19951) and iMehlert et al.l 
(2000J. Here we estimate the errors to be Alogrgff = 
0.1. This is slightly higher t han th e uncertainties given in 
lJ0rgensen. Franx & Kiaerg^dl (|1995|). b ut accounts for possi- 
ble svstematic errors (ISaglia et al.lll997l) . Stellar masses were 
computed from our best-fit ste llar mass-to-l ight ratios T and 
/?-band luminosities of Mehler fet al.l (120001) . In case the best 
fit is obtained with a lo garithmic halo, T is taken from column 
4 of Tab. 2 in IThomas et al. (2007a). In case of an NFW fit, 
T comes from column 8 of the same table. The best-fit stellar 
mass-to-light ratios of GMP3414 and GMP4822 are given in 
column 1 of Tab.|4](cf. App.|A]i. 

3. DARK MATTER SCALING RELATIONS 

Fig. [T] shows the scalings of dark halo core-radii r/, and 
halo asymptotic circular velocities v/, with B-band luminosity 
Lb and stellar mass (the corresponding galaxy parameters 
with errors are listed in Tab.[T]). Both, halo core-sizes and halo 
circular velocities tend to increase with luminosity and mass. 
The case for a correlation between v/, and Lb is weak, if the 
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Fig. 1. — Halo core-radius r;, and circular velocity v;, versus B-band luminosity Lb (a,c) and versus stellar mass M* 
(b,d). Large symbols: Coma ellipticals (red: central stellar population age tq > 6Gyr, orange: tq < 6Gyr, details in the 
text); thick solid lines: fits to galaxies with rp > 6Gyr; small symbols: round early-types from lGerhard et aU t20Q l): dot- 
ted: spiral galaxy sca ling relations from lPersic. Salucci & Stel ( 1996a-b); short-dashed: spiral galaxy scaling relations from 
IKormendv & Freeman (2004i) . 



sample as a whole is considered (cf. column 8 of Tab. |2]i. 
However, four galaxies (GMP0144, GMP0756, GMP1176 
and GMP5975) separate from the rest of the sample galax- 
ies in having both noticeably smaller r/, and v/,. These galax- 
ies are shown in light color in Fig. [T] As a general trend, 
halo parameters tend to scale more tightly with luminosity 
and mass when these galaxies are omitted. The solid lines in 
Fig. [T]show corresponding log-linear fits^. For comparison, 
in Tab. |2] we give both, fits to all Coma galaxies as well as 
fits to the subsample without the four galaxies offset in Fig.[T] 
The difference between these four galaxies and the rest of the 
sample is further discussed below. 

The logarithmic halos of equation ([U have two free param- 
eters. Any and pair of r/,, v/, or pi, characterizes a specific 
halo. Fig.|2]shows a plot of pi, versus r/,. Both halo parame- 
ters are clearly correlated. A linear relation fits the points with 
a minimum Xred - 0-41 (per degree of freedom; cf. Tab. |2]i. 
This rather low value partly derives from a degeneracy be- 
tween the halo parameters in the dyna mical modeling (e.g. 
iGerhard et al. 1998, Thomas et al. 2004) which correlates the 
errors in both quantities. In Tab.|2] such a correlation between 
the errors is not taken into account and the x%d might be thus 



^ F its for this paper are performed with the routine fitexy of IPress et alj 
(1991 . 



underestimated. A Xred much larger than unity would indicate 
some intrinsic scatter in Fig. |2] whereas the low Xred quoted 
in Tab. |2] formally rules out any intrinsic scatter. Note that 
dark matter halos in cosmological A^-body simulations can 
be approximated by a two-parameter family of halo models, 
where the parameters are corr elated qualitatively in a similar 
way as revealed by Fig.l^lCe.g. iNavarro. Frenk & Whitel[T996l 
IWechsler et al.ll2002l) . but with some intrinsic scatter 

The four galaxies offset in Fig. [T] are also slightly off- 
set in Fig. |2] However, given the large uncertainties, this 
is not significant and Fig. |2] is consistent with the halos of 
the four galaxies belonging to the same one-parameter fam- 
ily as established by the remaining Coma galaxies. This im- 
plies that the four galaxies primarily differ in the amount of 
stellar light (and stellar mass, respectively) that is associated 
to a given halo. Noteworthy, the fo ur galaxies have stellar 
ages To < 6Gyr (iMehlert et al.ll2003l) . while all other Coma 
galaxies are significantly older (mostly tq > lOGyr). A mere 
stellar population effect, however, is unlikely to explain the 
offset of the four galaxies. In this case differences to other 
Coma galaxies should vanish when galaxies are compared at 
the same stellar mass, which is not consistent with Fig.[TJ)/d. 
It should be noted, though, that the stellar masses used here 
are taken from the dynamical models. A detailed compari- 
son with mass-to-light ratios from stellar population synthe- 
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log Tj^/kpc 

Fig. 2. — Central dark matter density pi, versus halo core-radius ri,. Sym- 
bols and lines as in Fig.[T] 

sis models is planed for a future publication (Thomas et al., in 
preparation). 

In Fig.[3t we plot r/, against rgff . Larger core-radii are found 
in more extended galaxies. Three of the four galaxies with 
young cores are again offset. These three galaxies have sim- 
ilar Teff than other Coma galaxies of the same luminosity (cf. 
Fig- El')- This makes a higher baryon concentration unlikely 
to be the cause of their small halo core-radii. An exceptional 
case is GMP0756: it has a small core-radius, a ratio rijr^n 
which is typical for Coma galaxies with old stellar popula- 
tions and a relatively small reff. The scalings of r^a and r/, 
with luminosity in Coma galaxies with old stellar populations 
imply a roughly constant ratio r/,/reff « 3. Moreover, disk 
galaxies of the same luminosity show a similar ratio. 

The fact that all four galaxies offset in Figs. [1] and [3] ap- 
pear at projected cluster-centric distances D > 1 Mpc and have 
young stellar populations suggests that they may have entered 
the Coma cluster only recently. We will come back to this 
point in the next Sec. 13. 11 

3.1. Comparison with round and non-rotating early-types 

iKronawitter et all (120001) and iGerhard et all (120011) stud- 
ied dark matter halos of 21 nearly round (E0-E2) and non- 
rotating galaxies (K2000 in the following). K2000 galax- 
ies have similar luminosities Mb and half-light radii r^ff as 
ours, but the K2000 sample contains a mixture of field el- 
lipticals and galaxies in the Virgo and Forna x clusters. The 
dynamical models of Kronawi tter et al] (12000) differ in some 
respects from the ones described in Sec.|2]and this will be fur- 
ther discussed below. H owever, in their mass decomposition 
IKronawitter et alj (l2000h assumed the same halo profile as in 
equation ([l]), such that we can directly compare their halo pa- 
rameters to ours (cf. small black dots in Figs.[T]-|3]l. 

We find halo parameters of both samples in the same range, 
but Coma galaxies of the same Lb have on average larger halo 
core-radii than K2000 galaxies. However, the halos of the 
K2000 galaxies themselves are not different from the ones 
around Coma early-types, as both belong to the same one- 
parameter family (cf. Fig. |2]i. The main difference is that 
K2000 galaxies are brighter (and have higher stellar mass) 
than Coma early-types with a similar halo. Can this be an 
artifact related to differences in the dynamical models? 



iafl 

O 







11 



10.5 - 



^) 1 

• 

~ 1 1 1 1 1 1 1 1 1 




_ 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 I 

^ b) -.-t ' J 

• 

L > 1 

-, , 1 , 1 , 1 , , 1 1 , , , 1 , , 1 , 1 , 1 , , 1- 



-0.5 0.5 1 1.5 2 

log r ../kpc 



Fig . 3 . — Halo core-radius r/, (a) and 6-band luminosity (b) versus effective 
radius r^ff. La rge symbols: Coma galaxies; small dots: round galaxies from 
IGerhard et alj ^200 1ft . Lines in panel (a) as in Fig.n 

Many of the K2000 models are based on B-band photom- 
etry, while we used Rc-hwd images for the Coma galaxies. 
Elliptical galaxies become bluer towards the outer parts and 
Z?-band light profiles are slightly shallower than 7?-band pro- 
files. Likewise, mass profiles of galaxies are generally shal- 
lower than their light profiles, such that there might be less 
nee d for dark matter in B-band models than in /?-band mod- 
els. IKronawitter et al.l (12000) checked for this by modelling 
one galaxy (NGC3379) in both bands and found comparable 
results. The photometric data is therefore unlikely to cause 
the differences between the two samples. 

K2000 galaxies were modelled assuming spherical sym- 
metry. Not all apparently round galaxies need to be intrin- 
sically spherical. Neglecting the flattening along the line- 
of-sig ht can result in an u nderestimation of a galaxy's mass 
(e.g. iThomas et al.ll2007bft . Based on the average intrinsic 
flatte ning of ellipticals in the luminosity interval of interest 
here, ' Kronawitter et al.l (l2000l) estimated that the assumption 
of spherical symmetry should affect mass-to-light ratios only 
at the 10 percent level. We expect the effect on v/, to be corre- 
spondingly small. In addition, it is not obvious why spherical 
symmetry should enforce systematically small r^, such that 
the different symmetry assumptions are also unlikely to ex- 
plain the more extended cores and higher circular velocities 
in Coma galaxy halos. 

Since the shape of a galaxy is related to its evolutionary 
history, the round and non-rotating K2000 galaxies could 
be intrinsically different from the mostly flattened and ro- 
tating Coma galaxies. Structural differences could also be 
related to the fact that K2000 galaxies are located in a va- 
riety of environments, with less galaxies in high density re- 
gions like Coma. For example, stellar population models in- 
dicate that field ellipticals are on average younger and have 
more extended star fo rmation histories than cluster galaxies 
(iThomas et al.ll20()5a^ . But a mere difference in the stellar 
populations can not explain the difference between K2000 
and Coma galaxies, as in this case the scalings with M, should 
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be similar in both samples. This is ruled out by Fig.[TJ)/d. 

Noteworthy, most of the K2000 galaxies in Figs. [T]-|3] ap- 
pear similar to the four Coma galaxies with distinctly small 
r/, and v/,. As it has been discussed above, these galaxies may 
have entered the Coma cluster only recently and - in this re- 
spect - are more representative for a field galaxy population 
rather than being genuine old cluster galaxies. A consistent 
explanation for both, the offset between young and old Coma 
galaxies on the one side and the difference between old Coma 
galaxies and K2000 galaxies on the other would then be that 
field galaxies have lower r/, and lower v/, than cluster galaxies 
of the same stellar mass. Because there are more field galax- 
ies in the K2000 sample than in Coma, Coma galaxy halos 
would be expected to have on average larger cores and to be 
more massive (consistent with Fig. [Til. A larger comparison 
sample of field elliptical halos is required to conclude finally 
upon this point. 

3.2. Comparison with spiral galaxies 

Two independent derivations of dark matter scaling rela- 
tions for spiral galaxies are included in Figs.[T]-[3]through the 
dotted and dashed lines. Dotted line s show scaling relations 
from Persic. Salucci & Stell (11996a !'?). They are based on 
maximum-disk rotation-cur ve decompositions wit h the ha lo 
density from equation ([T]l. iPersic. Salucci & Stell ([T996allbl) 
give halo core-radii scaled by the optical disk radius. To 
reconstruc t the underlying rel ationship between r/, and Lb, 
we follow lGerhard et al] ( 120011) and assume exponential disks 
with 

!^^=8.4f-^V'' (4) 
^kpc^ \WLq) 



(the empirical fit in ' Gerhard et al.ll2001l has been transformed 
to our distance scale). 

Dashed lines in Figs. [T]-[3]fit the combined sample of 55 ro - 
tation curve decompositions of lKormendv & FreemanI (120041) . 
These authors discuss various halo profiles, but we here only 
consider non-singular isothermal dark matter halos. Though 
these are most similar to equation ([TJ, isothermal cores and 
circular velocities are not exactly identical as in logarithmic 
halos. To account for the difference, we fitted a logarithmic 
ha lo to a non-singular i sothermal density profile (cf. Tab. 4. 1 
of lBinnev & TremainelfT98 7). The fit was restricted to the re- 
gion with kinematic data (typically inside two core-radii). We 
found that the logarithmic halo fit yields a 3 percent larger 
core-radius r/, and a 10 percent smaller circular velocity v/,. 
The central halo density is reproduced to 0.001 dex (not sur- 
prising given that the cores of the two profiles were matched). 
Thus, under the assumption that fits performed with the two 
profiles indeed match inside two core-radii, a correction of the 
derived halo parameters is no t needed. Scaling relations from 
iKormendv & FreemanI (|2004|) are shown in Figs.[T]-|2]without 
a ny correction. 

IPersic. Salucci & Stell (Il996allbl) and lKormendv & FreemanI 
(|2004) discuss the scaling of disk galaxy halos with B-band 
luminosity. In order to compare early-types and spirals also at 
the same stellar masses, we used 



= 2x 



lOiiL. 



0.33 



(5) 



for the stellar mass-to-light ratios of disks. Equation (|5]l is 
derived fro m the Tully -Fisher and stellar-mass TuUy-Fisher 
relations of lBell & De JongI ( l2001h . 



Both in luminosity and in stellar mass, spirals and ellipti- 
cals follow similar global trends. However, while old Coma 
early-types have halo core-radii of similar size as spirals with 
the same B-band luminosity, the asymptotic halo circular ve- 
locities are 2.4 times higher than in corresponding spirals. In 
contrast, early-types with young central stellar populations 
have about 4 times smaller core-radii than spirals, but simi- 
lar asymptotic halo velocities. When galaxies are compared 
at the same stellar mass, then differences between ellipticals 
and spirals become larger (60 percent smaller r/, and 1 .8 times 
higher v/, in old Coma early-types compared to spirals; 90 
percent smaller r/, and 20 percent smaller v/, in Coma galax- 
ies with young central stellar populations). In addition, the 
halos of early-types and spirals do not belong to the same 
one-parameter family (cf. Fig.|2]i. At a given /-/, dark matter 
densities in ellipticals are about 0.5 dex higher than in spirals. 

4. THE DARK MATTER DENSITY 

Fig. m shows scaling laws for dark matter densities. The 
central dark matter density p/, (cf. equation O of the logarith- 
mic halo fits is plotted in panels (a) and (b) versus luminosity 
and stellar mass. For panels (c) and (d), the best-fit dark mat- 
ter halo of each galaxy (being either logarithmic or NFW) is 
averaged within 2 rgff : 



3 MDM(2reff) 
4^ (2reff)3 



(6) 



Here, Mdm('") equals the cumulative dark mass inside a sphere 
with radius r. Average dark matter densities are quoted in 
column (2) of Tab.|3| 

The general trend for both densities is to decrease with in- 
creasing galaxy luminosity. Thereby the central densities ph 
scatter more than the averaged (/Odm)- For two reasons, (pdm) 
quantifies the actual dark matter density more robustly than 
Ph. Firstly, our estimate of the very central dark matter den- 
sity depends strongly on the assumed halo profile. Instead, 
averaged over 2 reft differences between logarithmic halo fits 
and fits with NFW-profiles are small compared to the statis- 
tical errors (averaged over the Coma sample NFW-fits yield 
0.1-0.2 dex higher (pdm) than fits with logarithmic halos). 

Secondly, the most significant differences between /?/, and 
(Pdm) occur in the four Coma galaxies with distinct halos dis- 
cussed in Sec. [3] These galaxies have young central stellar 
populations. If the bulk of stars in these galaxies is old, how- 
ever, then the related radial increase of the stellar mass-to- 
light ratio could contribute to their small r/, and large p/,. This, 
because in our models it is assumed that the stellar mass-to- 
light ratio is radially constant. By construction then, any in- 
crease in the mass-to-light ratio with radius (being either due 
to a stellar population gradient or due to dark matter) is at- 
tributed to the halo component. In galaxies with a significant 
increase of the stellar M/L with radius, the 'halo' component 
of the model thereby has to account for b oth, additional stellar 
and possible dark mass (lThomasll2006l) . Any contamination 
with stellar mass will be largest at small radii, where the in- 
crease in the stellar M/L dominates the shape of the mass pro- 
file. The averaging radius in equation ^ is therefore chosen 
as large as possible. The value of 2 reff is a compromise for the 
whole sample, because the kinematic data extend to 1 - 3 r^a 
and the averaging should not go much beyond the last data 
point. 

Compared to the Coma galaxies, the majority of K2000 
galaxies have larger p/,. After averaging inside 2 reft, the halo 
densities in both samples become comparable, however. In 
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Fig. 4. — Central halo density pi, and average dark matter density (puu) inside Ir^ff versus luminosity Lb (a,c) and stellar 
mass M, (b,d). Solid lines: fits including all Coma galaxies; dotted and dashed lines: spiral galaxy scaling relations as in 
Fig.E] 



this respect, the K2000 galaxies again resemble the four Coma 
galaxies with young stellar cores. 

Dotted and dashed lines in Fig. |4] show spiral galaxies. 
Their halo densities need not to be averaged before compari- 
son, because core sizes of spirals (in the considered luminos- 
ity interval) are larger than Ir^ff. Averaged over the whole 
sample, we find dark matter densities in Coma early-types a 
factor of 6.8 higher than in spirals of the same luminosity. If 
early-types are compared to spirals of the same stellar mass, 
then the overdensity amounts to a factor of 13.5. Does this 
imply that spirals and ellipticals of the same luminosity have 
formed in different dark matter halos? 

5. BARYONIC CONTRACTION 

Even if ellipticals and spirals would have formed in simi- 
lar halos, then the final dark matter densities after the actual 
galaxy formation process could be different, since the baryons 
in ellipticals and spirals are not distributed in the same way. 
This effect can be approximated as follows: assume that (a 
spherical) baryonic mass distribution M*(r) condenses slowly 
out of an original haloH-baryon distribution M,(r). The halo re- 
sponds adiabatically and contracts into the mass distribution 
MuM(f)- If the original particles move on circular orbits then 
r[MJr)+MDM(r)] = nMi(ri) (7 ) 
turns out to be an adiabatic invariant jBlumenthal et alj l986). 

In case of the Coma galaxies, and Mdm are known from 
the dynamical modeling and equation (|7]i can be solved for 



Mi^. It characterizes the original halo mass distribution before 
the actual galaxy formation. Had a disk with baryonic mass 
grown in this original halo - instead of an early-type - 
then the halo contraction would have been different such that 
in general Mpjyj t^'Mdm- The difference between M^p^ and 
Mdm actually determines how much dark matter densities of 
ellipticals and spirals would differ, if both had formed in the 
same original halos. To quantify this further, let's consider 
the spherically averag ed mass distributi on of a thin exponen- 
tial disk for M° (e.g. iBlumenthal et al.l ll986). It is fully de- 
termined by a scale-radius and a mass. For a given Coma 
elliptical with luminosity Lg, the scale-radius and the mass of 
a realistic disk with the same luminosity can be taken from 
equations ^ and (|5]l. Then, given and M, (the recon- 
tracted Coma galaxy halo), equation (|7]i can be solved for the 
baryon-contracted halo M^y^ around the comparison disk (cf. 
IBlumenthal et al.lll986h . Once M^yi is known, the average 

(Pdm)° follows directly. 

If ellipticals and spirals (of the same Lb) would have formed 
in the same halos, then 



(Pdm) 
(Pdm) 



(8) 



' Because of the large core-radii in some galaxies (cf . Fig. [TJ.d) it is not 
always possible to find Mj for logarithmic halos. Therefore, here we only 
consider the best-fit NFW-halo of each galaxy. 
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Fig. 5. — Estimated formation redsiiifts Zform of Coma galaxy lialos versus luminosity without baryon correction (a) and 
with baryon connection (b). Symbols and colors as in Fig. [J 



should fully account for the observed ratio of elliptical to spi- 
ral dark matter densities. However, averaged over the Coma 
sample we find Shar ~ 2 and, thus, that the higher baryon con- 
centration in early-types is not sufficient to explain the factor 
of 6.8 between the dark matter densities of ellipticals and spi- 
rals at constant luminosity. 

In general, the observed dark matter density ratio (Jobs be- 
tween ellipticals and spirals will be a combination of a differ- 
ence in the halo densities before baryon infall and a factor that 
comes from the baryons. Let dha\o denote the baryon-corrected 
dark matter density ratio, then the simplest assumption is 

Sohs = 5bar X ^halo, (9) 

with (5bar from equation (O. After applying this approximate 
baryon correction, dark matter densities in Coma ellipticals 
are still a factor of 5haio = 3.4 higher than in spirals of the 
same luminosity. If the comparison is made at the same stel- 
lar mass, then (5haio = 6.4. Note that our baryonic contraction 
corrections are likely upper limits, because in equation (|5]l we 
only account for the stellar mass in the disk. In the presence 
of gas, the baryonic disk mass will be larger and so will be the 
halo contraction. The dark matter density contrast relative to 
the original elliptical will be therefore smaller. 

Concluding, the differences between the baryon distribu- 
tions of ellipticals and spirals are not sufficient to explain the 
overdensity of dark matter in ellipticals relative to spirals of 
the same luminosity or stellar mass. Ellipticals and spirals 
have not formed in the same halos. Instead, the higher dark 
matter density in elUpticals points to an earlier assembly red- 
shift. 

6. THE DARK-HALO ASSEMBLY EPOCH OF COMA EARLY-TYPES 

In order to evaluate the difference between elliptical and 
spiral galaxy assembly redshifts quantitatively, let's assume 
that dark matter densities scale with the mean density of the 
universe at the assembly epoch, i.e. pdm oc (l+Zform)^ (we will 



brm andzi,„ denote 



discuss this assumption in Sec. 17.2b . Let anu ^.^^^ 
the formation redshifts of ellipticals and spirals, respectively 
then 



1/3 



(10) 



(iGerhardet all 1200 Ih . with Pdm some measure of the dark 
matter density, e.g. Pdm = (pdm). Equation ( fTOl i can be solved 
for 



^form 



[1+4 



<5'/3-l. 



(11) 



where we have omitted the dependency of dark matter densi- 
ties and formation redshifts on Lb and defined 5 = p%yi/ Pom- 
Equation ( fTTT i allows to calculate formation redshifts of Coma 
and the observed 5. Two estimates based 



ellipticals from z|,^ 

on different assumptions about S and ^1,^^^^ are shown in Fig.|5] 



and are further discussed below. For each case, formation red- 
shifts were calculated with both disk halo scaling laws shown 
in Fig.|4]and the two results were averaged. 

Raw formation redshifts without any baryon correction (S = 
(Jobs) are shown in Fig.|5^ (and listed in column 3 of Tab. [3]). 
We considered a wide range of spiral galaxy formation red- 
shifts zfg^^ G [0.5,2] and the related uncertainty in zf^,^j„ is 
indicated by vertical bars. Our fiducial value is Zf„,.^ = 1, be- 
cause regular disks become rare beyond z > 1 dConselice et alj 
12005b . We did not allow for a luminosity dependence of spiral 
galaxy formation times. For the Coma ellipticals we then find 
4rm ranging from zf„,^ ~ 0.5 to z|^^ w 5, with the majority 
of galaxies having formed around zf^^^ ~ 3. Brighter galaxies 
have assembled later than fainter galaxies. 

Coma galaxy assembly redshifts shown in Fig.|5]3 (see also 
column 5 of Tab. [3]) include the baryon correction of Sec. |5] 
because we used S = (Jhaio in equation ( fTTT i. Moreover, be- 
cause fainter spirals have denser halos than brighter ones, we 
allowed for a luminosity dependent z|,rm(^)- ^^i analogy to 
equation ( fTOb assume 



l+4-.n(-L) ^f plM(L) \ 
1+4^(^-0) VPDM('Lo)y 



1/3 



(12) 



and z§ = 1 for a reference luminosity logLo/L© = 10.5 (the 
dashed line in Fig. [8]: illustrates the resulting zf^^.^ as a func- 
tion of L). Because (Jhaio ^ <^obs, the baryon corrected 
in Fig. |5l) are lower than the uncorrected ones in Fig. |5ji. 
The typical assembly redshift reduces to z^^^ ~ 2, as com- 
pared to z|,f„ ~ 3 without the correction. The baryon cor- 
rection is mostly smaller than the uncertainty related to our 
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Fig. 6. — As Fig. [5] but estimated formation redsliifts Zform of Coma galaxy halos are plotted versus the formation epoch 
z(tq) of the central stellar population (from column 7 of Tab.|3); dotted lines: one-to-one relations. Symbols and colors as in 
Fig.E 



ignorance about z^^^.^ (vertical bars in Fig. [Sji correspond to 
Zq G [0.5,2]). The trend for lower zf^j.^^ in brighter galaxies 
is slightly diminished by the baryon correction such that the 
dependency of on L in Fig.jS]? mainly reflects the lumi- 
nosity dependence of spiral galaxy assembly redshifts. 

Fig. |6] compares halo asse mbly redshifts with c entral stel- 
lar population ages tq from iMehlert et al.l (l2003l) . (We use 
Ho = 70kms"'Mpc"', D,\ = 0.75 and Vim = 0.25 to transform 
ages into redshifts.) Largely independent from applying the 
baryon correction or not, the agreement between the two red- 
shifts is fairly good for about half of our sample. Among the 
remaining galaxies, some have halos which appear younger 
than their central stellar populations. This could indicate that 
the stellar ages are overestimated (they are sometimes larger 
than the age of the universe in the adopted cosmology; cf. 
Tab. O. It could also point at these galaxies having grown 
by dry merging. In a dry merger, the dark matter density can 
drop, but the stellar ages stay constant. In Coma galaxies with 
young stellar cores, the halo assembly redshifts are instead 
larger than the central stellar ages. This indicates some sec- 
ondary star-formation after the main epoch of halo assembly. 

7. COMPARISON WITH SEMI- ANALYTIC GALAXY FORMATION 
MODELS 

In the following we will compare our results to semi- 
analytic galaxy formation models. To this end we have con- 
structed a comparison sa mple of synthetic ellipticals and spi- 
rals using the models of iDe Lucia & Bl aizol) (l2 0O7|). whic h 
are based on the Millennium simulation (Springel et al.l2005h . 
Comparison ellipticals are selected to rest in dark matter clus- 
ter structures with virial masses l arger than Myj > lO'^Mf:^ 
and to obey Mb bulge - Mb < 0.4 ( Simien & De VaucouleursI 
\l9m . We ignore galaxies at the centers of simulated clusters 
since we have omitted the two central Coma galaxies from the 
analysis in this paper Likewise, we exclude from the com- 
parison galaxies that have been stripped-off their entire halo, 
because the only Coma galaxy that possibly lacks dark matter 
inside 3reff has been excluded from the analysis in this pa- 
per as well (cf. Sec. |2|i. Isolated field spirals are drawn from 



objects with bulge > 1-56 in the semi-analytic models 
(ISimien & De Vaucouleurslll986h . 

Simul ated galaxies were chosen randomly from the cata- 
logue of De Luc ia & Blaizoi (|2007|) in a way such that each 
of six luminosity intervals (between Mb = -17 and Mb = -23; 
width AMb = 1 .0) contains roughly 50 galaxies. We use dust- 
corrected luminosities Mb of the semi-analytic models. 

7.1. Dark matter density 

Dark matter halos of simulated galaxies are reconstructed 
from tabulated virial velocities Vyir, virial radii ryir, and maxi- 
mum circular velocities Vj^ax as follows. It is assumed that the 
halos can be approximated by an NFW-profile (cf. equation 
[3]), in which case the circular velocity profile reads 

^ 1 ln(l + cx)-cx/(l+cx) 

' (13) 



Vcirc('") 



X ln(l+c)-c/(l+c) 

Here x = rjr^i^ and the halo concentration is defined by c = 
Tvii/ rj. The maximum circular velocity Vmax of a n NFW halo 
occurs at r « 2rvii /c jNavarro. Frenk & White! [19961) . such 
that (with equation[T3]) 



4.63 



ln(l + c)-c/(l+c)' 



(14) 



Using the tabulated Vvir and v^ax this equation can be numer- 
ically solved for the halo concentration c, which in turn de- 
termines rs = ?"vir/c and, thus, the entire NFW profile of the 
halo. 

Before compared to the Coma galaxy models, halo densities 
are averaged within 2 rgff (cf. equation|6l). In case of simulated 
spirals we use effective radii from the empirical relation (|4|. 
For ellipticals, we assume 



kpc 



= 15.34 



1.02 



(15) 



which is a fit to the Coma data. 

Fig.|2]shows that the average dark matter densities {puu) of 
the Coma early-types match fairly well with semi-analytical 
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Fig. 7. — Average dark matter density (pdm) versus Lb in simulated cluster ellipticals (a ) and in simulated field spir als (b). 
Large symbols and lines as in Fig.[T] Simulated galaxies from the semi-analytic models of lDe Lucia & Blaizoll 120071) . 



models. This is remarkable, because the simulations do not 
take into account the halo response during baryon infall. 
Therefore, either the net effect of the baryons on the dark mat- 
ter distribution is small in the analyzed population of galax- 
ies or there is actually a mismatch between the halos of ob- 
served galaxies and the A^-body models. It may also be that 
real galaxies do not have maximum stellar masses. This can 
be checked by the comparison of dynamically derived stellar 
mass-to-light ratios with independent stellar population syn- 
thesis models (Thomas et al., in preparation). 

Similarly to what is found in real galaxies, the dark mat- 
ter densities of spirals are lower than in ellipticals in semi- 
analytic models (cf. Fig. |7]), but the density contrast in ob- 
served galaxies is larger. Again, a major uncertainty here is 
that the simulations do not take into account the gravitational 
effect of the baryons. 

7.2. Assembly redshift 

Formation redshifts of simulated and observed galaxies are 
compared in Fig. [8] Coma galaxy Zfonn are from Sec. |6] and 
both cases discussed there - with and without baryon correc- 
tion - are shown separately in panels (a) and (b), respectively. 
Formation redshifts of simulated galaxies are defined as the 
earliest redshift, when a halo has assembled 50 percent of its 
mass. Since we are mainly interested in cluster ellipticals, 
we need to take into account that interactions between the 
cluster halo and a galaxy's subhalo cause a mass-loss in the 
latter Although cluster-galaxy interactions happen in both 
simulated and observed galaxies, the mass-loss in the simu- 
lations may be overestimated because of the finite numerical 
resolution and the neglect of the baryon potential. In partic- 
ular, for simulated subhalos with very low masses at z = 
the derived formation redshifts may be artificially high, when 
defined according to the assembly of half of the final mass. 
To avoid such artificially large assembly redshifts, we define 
Zfoim of simulated galaxies as the earliest time when half of the 
maximum mass was assembled, that a single progenitor in the 
merger tree of given galaxy had at some redshift. Our assump- 
tion is that even if dynamical interactions between cluster and 
galaxy halos take place, they do not significantly affect the 



very inner regions < 2 r^a of interest here. In case of field spi- 
rals, formation redshifts defined either from the final or from 
the maximum mass are very similar 

Without a baryon correction, our estimates of Coma galaxy 
formation redshifts are on average higher than in the semi- 
analytic models (Fig. [8^). This, although (1) the dark mat- 
ter densities of ellipticals match with the simulations and 
(2) our assumption about the formation redshifts of spirals 
(^form ~ 1) is consistent with the simulations. The origin for 
the offset between Coma galaxies and semi-analytic models 
in Fig. [8^ is that the density contrast between halos of el- 
lipticals and spirals is larger in observed galaxies than in the 
simulations. After applying the baryon correction, the Coma 
galaxy formation redshifts become consistent with the simu- 
lations (Fig.lSJ)). This result indicates that the discrepancy be- 
tween the measured and the simulated density ratio Pom/Pdm 
is due to baryon effects. 

Our Coma galaxy formation redshifts are based on the as- 
sumption that (pdm) (1+Zform)"'- Fig.|9]shows (pdm) versus 
(1 +Zform) explicitly. Independent of including a baryon cor- 
rection or not, the slope of the relationship between (pdm) and 
( 1 +Zform) in the Coma galaxies is roughly parallel to simulated 
A^-body halos. This confirms that our assumption for the scal- 
ing between (pdm) and Zfon-n is approximately consistent with 
the cosmological simulations. 

Concerning the absolute values of the dark matter densities 
it has been already stated above that they are only consistent 
with the simulations if either the net effect of the baryons is 
zero in the case of ellipticals or if galaxies do not have max- 
imum stellar masses. The former case would imply that ha- 
los of spiral galaxies experience a net expansion during the 
bar yon infall (several processes h ave been proposed for this, 
e.g. Binnev, Gerhard & Sil^lldoJ). 

De Lucia et al. (2006) quote a stellar assembly redshift be- 
low z < 1 for simulated ellipticals more massive than > 
1O"M0. The halo assembly redshifts in Fig. [8] are mostly 
above z > 1. In part, this is due to the fact that we only 
consider semi-analytic galaxies in high-density environments 
similar to Coma. In addition, formation redshifts defined ac- 
cording to the stellar mass assembly and the halo assembly. 
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Fig. 8. — Comparison between assembly redshifts Zfonn of simulated and observed galaxies (symbols as in Fig.|7). Both formation redshift 
estimates for Coma galaxies with and without baryon correction are shown (pan els a ,b). For comparison, panel (c) shows formation redshifts 
of sim ulated spirals (dots) and of observed spirals (dashed-line; cf. equation 112) . Simulated galaxies from the semi-analytic models of 
iDeT ucia & Blaizot (2007). 



respectively, are not always equal. For example, in our com- 
parison sample of simulated cluster ellipticals we find an av- 
erage dark halo assembly redshift (zform) = 1-50 for galaxies 
more massive than > IO^Mq. Evaluating for the same 
galaxies the redshift Zf^j.^^ (when half the stellar mass is as- 
sembled) yields (Zfoi-n,) = 1.07. That Zhnn > Zform is plausible 
if some star formation is going on between < z < Zform in the 
progenitor and/or in the subunits that are to be accreted after 
Zform- It should also be noted that the simulations do not take 
into account stellar mass-loss due to tidal interactions. 

8. SUMMARY 

We have presented dark matter scaling relations derived 
from axisymmetric, orbit-based dynamical models of flat- 
tened and rotating as well as non-rotating Coma early-type 
galaxies. Dark matter halos in these galaxies follow similar 
trends with luminosity as in spirals. Thereby, the majority of 
Coma early-types - those with old stellar populations - have 
halo core-radii r/, similar as in spirals with the same B-band 
luminosity, but their asymptotic halo velocities are about 2.4 
times higher In contrast, four Coma early-types - with young 
central stellar populations - have halo velocities of the same 
order as in comparably bright spirals, but their core-radii are 
smaller by a factor of 4. Differences between spirals and el- 
lipticals increase, when the comparison is made at the same 
stellar mass. The average halo density inside Ir^ff exceeds 
that of comparably bright spirals by about a factor of 6.8. If 
the higher baryon concentration in ellipticals is taken into ac- 
count, the excess density reduces to about a factor of 3, but if 
ellipticals and spirals are compared at the same stellar mass, 
then it is again of the order of 6.5. 

Our measured dark matter densities match with a 



comparison sample of simulated cluster ellipticals con- 
structed from the semi-an alytic galaxy formation models of 
iDe Lucia & Blaizot (|2007[) . These synthetic ellipticals have 
Zform « 0.5 -4 and higher dark matter densities than simulated 
field spirals, which are appear on average around z^j.^^ ~ 1 . 

Assuming for local spirals Zf^^^^ = 1 as well, and assuming 
further that the inner dark matter density scales with the for- 
mation redshift like (1 +Zfomi)"^, our results imply that ellip- 
ticals have formed Azform ~ 1 - 2 earlier than spirals. With- 
out baryon correction, we find an average formation redshift 
around ztrnm ~ 3, which is slightly larger than in semi-analytic 
galaxy formation models. Accounting for the more concen- 
trated baryons in ellipticals, the average formation redshift 
drops to Zfoi-m ~ 2. 

For about half of our sample, dark halo formation red- 
shifts match with co nstraints derived from stellar populations 
dMehlert et al.l2003l) : the assembly epoch of these (old) early- 
types coincides with the epoch of formation of their stellar 
components. 
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APPENDIX 

GMP3414 AND GMP4822 

The best-fi t model parameters for the galaxies GMP3414 and GMP482 2 (which were not inc luded in the original sample of 
iThomas et al.. 2007a) are given in Tab. |4] The table is similar to Tab. 2 of iThomas et al.l (l2007al) and we refer the reader to this 
paper, in case more detailed information about the parameter definitions are required. The best-fit models with and without dark 
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Fig. 9. — Average halo density (pdm) versus assembly redshift. Large symbols: Coma galaxies without baryon correction 
(a) and with bary on conection (b). Small sy mbols: simulated cluster ellipticals and simulated field spirals from the semi- 
analytic models of lDe Lucia & Blaizoll 120071) . 

matter halo are compared to the observations in Figs. [TOl and [TTl In both galaxies the best-fit inclination is / = 90°, but the 68 
percent confidence regions include models at / > 70° (GMP3414) and / > 50° (GMP4822). 

GMP3414 and GMP4822 were observed with the Wide Field Planetary Camera 2 (WFPC2) on board the HST as part of the 
//5r proposal 10844 (PI: G. Wegner). For each galaxy two exposures with 300s each were taken with the filter F622W. Four 
other objects were previously ob served as part of this proposal and a full description of the respective observational parameters 
and the data analysis is given in ICorsini et alJ (l2008l) . A table with the final photometric parameters will be published in the 
journal version of this paper. 
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GMP3414 




GMP3414 




Fig. 10. — Upper panel: joint ground-based and HST photometry of 
GMP3414/NGC4871. Lines: best-fit deprojection (red) and its edge-on reprojection 
(blue). Lower panel: stellar kinematics along major axis (left/red), along the minor axis 
(middle/blue) and along a third axis parallel to the major axis with an offset of reff/2 
(right/green); filled and open circles refer to different sides of the slits; dotted: best-fit 
model without dark matter. 
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Fig. 1 1 .— As Fig.[T0]but for GMP4822/NGC4841A. 
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TABLE 2 
Dark matter scaling relations. 
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fits omitting galaxies with young stellar cores 



kpc 

Ih_ 
kpc 

v> 

km s^' 

I'fi 
km s^' 

Ph 

Ih_ 
kpc 

Ph 
Mqpc^ 
Ph 

(PDM> 
Mq pc3 

(pdm) 



1.54±0.21 


0.63 ±0.33 


0.16 


0.16 


0.44 


0.002 




0.98 ±0.14 


0.54 ±0.29 


0.18 


0.17 


0.44 


0.007 




0.78 ±0.11 


0.21 ±0.21 


0.14 


0.07 


0.24 


0.016 




0.59±0.10 


0.19±0.17 


0.12 


0.07 


0.24 


0.008 


m 


0.68 ±1.35 


-1.64 ±0.96 


0.04 


0.18 


0.75 


0.006 


m 


0.79 ±0.21 


0.62 ±0.30 


0.10 


0.14 


0.44 


0.017 


m 


-1.81 ±0.36 


-1.02 ±0.54 


0.31 


0.34 


0.75 


0.028 




-0.94 ±0.17 


-0.81 ±0.43 


0.31 


0.36 


0.75 


0.090 




-2.01 ±0.20 


-1.04 ±0.30 


0.77 


0.31 


0.41 


0.012 




-1.12±0.11 


-0.74 ±0.25 


0.86 


0.33 


0.41 


0.020 





Note. — (1,2) Fitted quantities; (3,4) parameters of linear fit with errors; (5) reduced Xred 
of the fit; (6) rms-scatter in logy; (7) mean error {/S.\ogy)\ (8) significance of the relation 
(probability V that there is no relation according to a Spearman rank order correlation test); (9) 
figure in which the relation is shown. 
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TABLE 3 

Halo densities, halo assembly redshfits and stellar ages 
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Note. — (1) Galaxy id from iGodwin, Metcalfe & PeacH ( 119831) : (2) average dark matter 
density (pdm) inside Ir^ff; (3) halo assembly redshift Zform according to equation dl It with 
5 = Sobs and Zfo,.,^ = 1; (4) as column (2), but including the baryonic correction Jbai defined in 
equation ([8} (we do not derive an error estimate for the baryonic contraction, but use the same 
errors in columns (2) and (4), respectively); (5) Zfonn as in column (3), but with baryon corrected 
S = (?ha io and Zforra(^) from equation ( 112b : (6) central stellar age to from Tab. B. 1 of Mehlert e t all 
( 12003 ^ (GMP3958 has no age estimate because of its very low H/3); (7) formation redshift z(ro) 
of the stars derived from column (6). In some galaxies the stellar age or its upper limit exceed 
the age of the universe in the adopted cosmology. In such cases only a lower limit is given for 
z(to) or the upper redshift error is set equal to 00, respectively. 



TABLE 4 

Model parameters for GMP3414 and GMP4822. 
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Note. — (I) galaxy id; (2) type of fit (SC: without 
dark matter; LOG: logarithmic h alo; NFW: halo profile 
from lNavarro. Frenk & W hite" 1996): (3) best-fit stellar T 
[Mq/Lq] (i?c-band); (4,5) best-fit logarithmic halo pa- 
rameters ri, [kpc] and vi, [kms~'l; (6,7) best-fit NFW c on- 
centration c and flattening q (cf. iThomas et al ] |2007al for 
details); (8) achieved goodness-of-fit (per data point). 



